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Abstract. This paper presents a comparative 
overview of the disciplines of modal testing 
used In structural engineering and system 
Identification used In control theory. A list 
of representative references from both areas Is 
given, and the basic methods are described 
briefly. Recent progress on the Interaction of 
modal testing and control disciplines Is dis- 
cussed. It Is concluded that combined efforts 
of researchers In both dlsclpl Ines are required 
for unification of modal testing and system 
Identification methods for control of flexible 
structures. 


Modal testing In the field of structures means 
the process of measuring signals produced by a 
structure and Identifying modal parameters 
(damping, frequencies, mode shapes and modal 
participation factors). System Identification 
In the field of controls means the process of 
measuring signals produced by a system and 
building a model to represent the system for 
control design. Techniques to Identify a model 
from measured data typical ly contain two steps. 
First, a family of candidate models Is chosen 
and then the particular member In this family 
Is determined wh tch satisfactorily describes 
the observed data based on some error crite- 
rion. If the Identified model Is a linear 
model In state space representation, the el gen- 
sol utlon of the model provides eigenvalues and 
eigenvectors that, in turn, determine modal 
parameters for structures. Correlation between 
the fields of modal testing and system Identi- 
fication for controls Is evident. 


The area of modal testing Is a well-developed 
dlsclpl Ine with strong exper 1 mental foundations 
[1-76]. The area of system Identification for 
control s I s we I I -dev el oped w I th sol Id theoretl - 
cal and methodological foundations [77-170]. 
While the development of each Individual area 
continues, there Is a need to provide a compre- 
hensive yet coherent unification of the areas. 
Active control of flexible structures will 
require the combined efforts of researchers In 
both disciplines. />mong these challenges Is 
control of large space antennas and platforms. 

* NASA Langley Research Center, Hampton, Virginia 


The areas of modal testing and system Identifi- 
cation encompass a multitude of approaches, 
perspectives and techniques whose Interrela- 
tionships and relative merit are difficult to 
sort out. As a result, It Is difficult for a 
nonspecl al I st to extract the fundamental con- 
cepts. It may take considerable effort to gain 
enough Intuition about a particular technique 
~o be able to use It effectively In practice. 

"he objective of this paper Is to present an 
overview of the paral lei historical development 
of modal testing used In structural engineering 
ond system Identification used In control 
"heory. A list of principal references is 
provided for studying the similarities and 
differences among the many approaches in both 
areas. 


MODAL TESTING 


This section contains a synopsis of the field 
of modal testing. The following three Items 
are providedi (1) a concise, yet complete, 
chronology of key developments that have oc- 
curred over the 40-year history of modal test- 
ing, (2) a brief summary of currently used 
approaches, and (3) a chronological reference 
list of key publications. 

The subject of modal testing has evolved con- 
tinuously since the 1940* s, and an extensive 
literature has been generated (for example, 
note the bibliographies on pp 1659-1734 In 

Proceedings of the 4.th International — Moda l 

Analysis Conference . February 1986). It Is 
beyond the scope of this paper to discuss all 
this activity. Interested readers are referred 
to these references, or to several other more- 
complete overviews which ha/e been written 
recently [51,53,64,74,75], for additional 
Information. 


Chronology of Key Developments. The chronology 
Is divided Into three separate eras: pre-1970, 

1970-1979 and 1980-present. The start of the 
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second and third eras accompanied significant 
Improvements In computer technology. The 
second era began with the widespread Introduc- 
tion of minicomputers (e.g., PDP-11) and the 
third with the avail ab! I Ity of much larger and 
faster computers for laboratory data analysis 
(e.g., VAX or mainframes). 

In the first era of modal testing, from ap- 
proximately 1940 through 1969, analog tech- 
niques were used almost exclusively. 
Laboratory computers were noT yet available, 
and mainframe computers were used only rarely 
for data analysis. Much of the earliest work 
occurred In the aircraft Industry, where test- 
ing was conducted to check the accuracy of 
calculated normal modes used In flutter and 
dynamic loads predictions. The two most sig- 
nificant contributions from this period were 
the works of Kennedy and Pancu Cl], who Intro- 
duced "circle fitting 1 for decomposing fre- 
quency response functions (FRF f s) Into the 
constituent modes, and Lewis and Wrlsley [2], 
who described a systematic approach for tuning 
Individual modes using multiple shakers and 
apportioned sinusoidal excitation. These two 
techniques, with numerous variations, were used 
In the majority of modal tests conducted prior 
to 1970. Testing was very time constmtng, 
however, and required considerable practice and 
skill for success. All I aboratory equipment 
was analog, and most data analysis was per- 
formed by hand. 


The second era occurred during the 1970 f s f 
sparked by the Introduction of laboratory 
minicomputer systems [15, 19] and the fast 
Fourier transform algorithm to compute fre- 
quency response functions [18, 30], Compared 
with the classical method of slowly sweeping a 
sinusoidal signal to generate FRF f s, these 
systems offered tremendous speed advantages. 
They were widely adopted by the modal testing 
community, with the exception of those organi- 
zations that had already made large Investments 
In multiple-shaker sinusoidal testing equip- 
ment. The minicomputer In these laboratory 
systems was used not only to compute FRF*s, but 
also for curve fitting the FRF data to estimate 
modal parameters. Many of the analysis tech- 
niques used during this period, however, were 
simply digital versions of techniques developed 
earlier (e.g., circle fitting or phase separa- 
tion techniques). This situation began to 
change during the second half of the 1970 ! s in 
conjunction with the availability of more 
powerful processors [3 1 D . Also, the use of 
mainframe computers for data analysis was 
beginning to occur [34,35]. There was consid- 
erable hesitancy to use mainf runes, however, 
because of the expense and because most analy- 
sis techniques In use at the time required 
considerable Interaction with the user. 


The third era of modal testing began around 
1930, again In conjunction with improved com- 
puters and data acquisition equipment. One of 
the most significant changes occurred In the 
use of mul tl pi e- Input random excitation [42], 
rather than the single-input random approach 
previously used most often. Mul 1 1 pi e- I nput 
excitation provides several advantages: (1) it 

Is consistent with mul 1 1 pi e-ref erence (multi- 
ple-input) modal Identification algorithms 
[49,55,60,61,67] which allow closely spaced 
modes to be better Identified, (2) It minimizes 
troublesome shifts In frequency and mode shapes 
which can occur when exciters are moved to 
different points on a structure, and (3) It Is 
significantly faster whenever several different 
locations and directions of excitation are 
needed to excite all the modes of Interest, 
which Is usual ly the case. In conjunction with 
mil 1 1 pi e- 1 nput methods. Improved methods of 
excitation [54,56], signal processing [65], 
frequency response calculation [59], and quick 
data analysis [66,69] have also been Intro- 
duced. Renewed attention Is being given to 
ARMA-type analysis techniques [72] which were 
studied years earl ier [11,13,35] and general ly 
thought to be too computational ly extensive for 
the volume of data obtained In modal tests 
(often >20 modes and >100 measurements). Faster 
and more powerful computers now make these 
approaches more practical. Faster computers 
also permit nonl I near It! es to be better de- 
tected and Identified, using new signal proc- 
essing techniques [57]. Nonl i near Itl es are 
also being better quantified using modern 
slepped-slne excitation techniques. In conjunc- 
tion with new modal Identification algorithms 
which can use unequal ly spaced frequency data 
[70]. Renewed Interest Is also occurring in 
the classical forced normal mode approach to 
modal testing. However, today the excitation 
tuning process Is being computerized more than 
ever [76], 

Current Methodology. As discussed above, modal 
testing methodology has changed considerably 
over the years — and It continues to change. 
There are now many different ways to conduct 
testing and data analysis. Essentially all 
techniques work wel I with simple structures, 
yot significant differences occur when used on 
complex, built-up structures. It Is still 
often difficult to deduce the exact source of 
these differences, however, because the "true" 
answers are unknown with experimental data. 
Better methods for comparing various identifi- 
cation techniques with complex data are needed 
[M], as are methods for conducting more real- 
istic and thorough simulations of complex 
structural behavior to generate data for these 
studies. To a large degree, the techniques 
used most often today are those which ha/e 
demonstrated repeatedly their tolerance of the 
complexities of real data. 
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Figure 1 provides an estimate of the general 
types of excitation and analysis now being used 
In the modal testing community In the United 
States (the percentage of multi- Input tuned 
sine dwell testing Is much higher In Europe). 
Mul tl pi e- 1 nput random excitation Is new very 
popular. It Is estimated that 45$ of all 
laboratory modal tests are being conducted this 
way. Still very popular 1 s the si mpl er, tradi- 
tional single- 1 nput random approach (3$). The 
percentage of people using multi- Input tuned 
sine dwell (forced normal mode) testing has 
been fal ling since around 197 5 and Is estimated 
to be about 10%. Of these, approximately 
one-half are now using computers to generate 
optimum tuning patterns. A smal 1 percentage of 
tests (5$) are also conducted using natural 
ambient excitation forces (e.g., wind on build- 
ings, waves on of f shore o I ! platforms), usually 
because artificial excitation Is Impractical. 
In many cases, the Identification results from 
these tests are much more ambiguous than those 
from control I ed-e>c1 tatlon tests. 
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Figure 1. Current U.S. Modal Test Methodology 
Preferences (Estimated). 

In the vast majority of testing (80$), fre- 
quency response functions are computed prior to 
using a modal Identification algorithm. There 
are many reasons why FRF f s are stll I generated 
so often, although there are new many Identifi- 
cation algorithms available which can analyze 
free- or forced-response time histories di- 
rectly. A significant factor Is tradition: 
experienced modal testing personnel can deduce 
considerable Information simply by observing 
frequency response functions. If time hlsto- 
rles are processed directly by the modal Iden- 
tification algorithm, many traditional 
eval uatlon cr Iterl a are unaval I abl e. The I nter- 
medl ate step of cal cul atl ng FRF f s will probabl y 
be skipped more often In the future as more 
experience and confidence Is gained with the 
newer, direct analysis approaches. 


For modal Identification, mul t l-degree-of- 
freedom (MDOF) time-domain algorithms are being 
ised the majority of the time (40$), followed 
by MDOF f req uency-doma! n (30$), and faster 
j \ ng I e— degree- of— freedom (SD0F) algorithms 
1 20$). In most cases, the data used by the 
MDOF time-domain algorithms are Impulse re- 
sponse functions obtained by Inverse Fourier 
Transformation of FRF' s. Only In a few l se- 
ated Instances are f ree-response data being 
nsed directly. MDOF time-domain algorithms 
:,tll I appear to be unsurpassed In their abl I Ity 
to analyze wideband data with many modes, while 
"heir frequency- domain counterparts ha/e been 
observed by seme people to be more accurate 
over narrow frequency bands, or when damping Is 
relatively high (e.g., greater than 5$). 

SYSTEM IDENTIFICATION USED IN CONTROL THEORY 

i) Ince the mld-1960's the f lei d of system lden- 
Mflcatlon has been an Important discipline 
within the automatic control area. One reason 
s the requirement that mathematical models 
within a specified accuracy must be used to 
,: 5 p p | y modern control methods. Another reason 
is the aval I abl I Ity of digital computers that 
can perform complex computations. It Is not 
the aim of this paper to present a detailed 
analysis of achievements In this continuously 
growing field or to give the state-of-the-art. 
See references [77-110] for more details. 

r requency and time domain methods give comple- 
mentary views of many Important problems In 
I Inear system theory and control theory. Some- 
times, the two methods ha/e been seen as ri- 
vals, particularly on Issues of Implementa- 
tion and appl Icatlon to real systems. Histori- 
cal ly, frequency domain methods dominated 
theory and practice of system Identification In 
control engineering prior to the 1960's. 

Frequency Domain Approach. Frequency domain 
Identification In control engineering gained 
relevance with stabll Ity and design methods 
based on frequency response measurements. 
Frequency response estimation began with the 
technique known as transfer function analysis. 
The sinusoidal transfer function analyzer Is 
recognized as a robust and practically useful 
non- parametr Ic Identification method. This Is 
due to the Intrinsic rel fabfl Ity with which the 
sinusoidal transfer function analyzer Is able 
to reject low-frequency drift and harmonic 
distortion due to nonl inearlty. However, the 
sinusoidal transfer function analyzer requires 
long test times to sequentially Identify each 
relevant point on a frequency response curve. A 
more widely used class of techniques has been 
developed around digital spectral analysis and 
numerical Fourier transforms. The Impact of 
the Fast Fourier Transform (FFT) developed by 
Cooley and Tukey Cl 13] upon digital and analog 
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spectral analysis was enormous. The FFT actu- 
al ly encompasses a whole family of algorithms, 
many of which are Included In the Institute of 
Electrical and Electronics Engineers (IEEE) 
collected reprint written by Rablner and Rader 
[117]. Modern digital spectral analysis, real 
time or off-line. Is normally achieved by the 
direct method which takes the FFT of data 
blocks and then averages the resulting spectral 
estimates. 

The direct estimation scheme can also be ap- 
pl led to determine estimates of frequency 
response functions by using a closed loop 
system design. The normal open loop system 
design may give biased results with physical ly 
unrealizable frequency response estimates. A 
closed loop system design may Improve this 
situation, leading to gain and phase estimates 
which are well behaved, provided that careful 
attention Is paid to windowing and al lasing. A 
more precise discussion of the accuracy of 
closed loop estimates Is given In Daval I [1 19], 
and Wei I stead [1 22]. 

The maximum entropy method (^€M) of spectral 
analysis was originated by Burg Cl 14] for 
analyzing geophysical data, and further devel- 
oped by Abies [118] and Ulrych and Bishop 
[120]. The basic philosophy of the method Is 
to construct spectral estimates that are con- 
sistent with all relevant data and are maxi- 
mally non-committal with regard to unavailable 
data. The maximum entropy method general ly 
gives superior spectral resolution to tradi- 
tional methods at the expense of Increased 
variability, and possibly erroneous splitting 
of spectral lines observed by Fougere, etc. 

[12137 

Frequency domain Identification, which, In the 
past, emphasized non-parametr Ic Identification, 
I.e., frequency response estimation, has lost 
popularity In recent years. This Is due to the 
fact that current control synthesis and design 
tools require parametric system models such as 
a state space representation, stochastic dif- 
ference equation or general Ized regression 
model. If a spectral analysis or transfer 
function analysis experiment Is conducted, 
least squares can be used to f It a parametric 
frequency response model by assembl Ing a set of 
N-measured frequency response points and solv- 
ing the unknown coefficients of the transfer 
function model [112,115,123]. However, these 
fitting methods may produce poor models which 
may not represent the underlying system well 
enough for controller design, particularly for 
complex structural dynamics problems. To a 
limited degree, these shortcomings are being 
investigated by current work on combined Iden- 
tification and control. Indeed, the pole-zero 
assignment technique developed by Wei I stead, et 
al. [124] seems to offer qualities of robust- 


ness. The technique developed by Juang and 
Suzuki [168], which uses estimated frequency 
spectra to Identify a state space model in 
modal space via system real Ization theory, also 
seems to offer a good model for control design. 

Tine Domain Approach. The time domain approach 
has dominated the control engineering litera- 
ture on system Identification over the past 20 
years. In this section, a general description 
o~ commonly used linear system Identification 
methods will be given. Time domain approaches 
are categorized according to the choice of 
model and the choice of Identification crite- 
rion for evaluating the estimation qual Ity. 
Basic methods are given In the references. 

Tae origin of the least squares method can be 
traced to Gauss [125] who formulated the basic 
concept and used It practically for astronomi- 
cal computation. Since then. It has been 
widely appl led to many problems. The recursive 
algorithm to calculate the least squares esti- 
mate has apparently been found Independently by 
several authors. The original reference seems 
to be Plackett [128]. An early and thorough 
treatment of the least squares method appl led 
to dynamic system Identification Is given by 
Astrom [94]. The statistical background for 
stochastic approximation was developed by 
Robbins and Monro [129]. Stochastic approxima- 
tion methods have also been derived by Sakrlson 
[138]. Computational algorithms are based on 
stochastic gradient methods for I inear regres- 
sion models. In LJung [161] the stochastic 
approximation approach Is used to derive recur- 
sive Identification algorithms for problems 
ether than I Inear regression models. 

The Kal man-Bucy filter [132] Isa state estima- 
tor. It Is Mayne [134] who draws attention to 
extending the Kal man-Bucy filter for parameter 
estimation of a state space model. The basic 
Idea of the instrumental variable method (e.g., 
Kendal and Stuart [133] and Young [143]) Is the 
generation of an extra signal, I.e., the In- 
strumental variable, which Is correlated with 
the useful signals of the process but which Is 
mcorrel a ted with noise. This el Imlnates the 
Mas error associated with least squares esti- 
mation. Recursive Instrumental variable meth- 
ods have been used extensively by Young [153]. 

The characteristics of noise corrupting the 
output of the system may not be wel I knewn. In 
the generalized least squares method, the 
parameter estimates may include estimates of 
noise parameters. Inspired by Clarke f s algo- 
rithm [139] for generalized least squares 
analysis, a recursive method was suggested by 
Hast l ng- James and Sage [142]. The principle of 
extended least squares Is that the calculation 
of the error between the true output and the 
estimated output Is based on past estimates of 
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system and noise parameters. The extended 
least squares algorithm was Independently 
derived by Panuska [1 40] and Young L 1 41 3 and 
widely used and rediscovered by Talmon and van 
den Boom [149]. In many practical problems of 
parameter estimation the problems arises of 
solving an overdetermined Ill-conditioned set 
of algebraic equations. To circumvent this 
problem, the error covariance matrix can be 
propagated In a square root form so that the 
positive sem I-def I n I te nature of the error 
covariance Is maintained to minimize the com- 
plexity of the statistical properties of the 
error estimates. A survey of square-root 
filtering techniques was given by Kaminski, et 
al. [146], 

The basic Idea of maxi mum- I Ike I Ihood estimation 
Is to construct a function of the data and 
unknown parameters cal led the I Ikel Ihood func- 
tion. The likelihood function Is essentially 
the probability density of observations. The 
estimate of parameters Is then obtained as the 
parameter set which maximizes this function. 
The method of maximum I lkel Ihood was developed 
by Fisher [126,127] although the basic Idea 
dates back to Gauss [125]. In the Bayesian 
approach, the parameters themselves are treated 
as a random variable. Based on observations of 
other random var 1 abl es that are cor rel ated w Ith 
the parameter. Information about Its value can 
be Inferred. Therefore, the parameter estimate 
Is expressed In terms of the probability dis- 
tribution conditioned by past history. 

Ho [135] showed that the Instrumental variable 
method, generalized least squares and extended 
least squares are closely related to stochastic 
approximation and Kalman-Bucy filtering. Actu- 
al ly , square-root f 1! terlng belongs to the same 
family, as wel I, but the computational concept 
differs essentially. For linear systems and 
Gaussian noise, the maximum- I Ikel 1 hood approach 
yields the same conditions for the parameter 
calculation as the least-squares approach. 
Although the preceding techniques have been 
widely used In the field of controls, formal 
dl rect appl Icatlon to modal parameter Identifi- 
cation for flexible structures has been mini- 
mal. Two techniques which has been extended 
and applied for modal Identification of struc- 
tures — minimum realization and lattice fil- 
tering — are discussed In the following 
section. 

Techniques Related to Modal Testing. In the 

field of controls, the process of constructing 
a state space representation using experimental 
data Is cal led system real Izatlon. A minimum 
real Izatlon Is a model with the smal lest state 
space dimension among models real Ized that have 
the same Input-output rel atlons w Ith 1 n a speci- 
fied degree of accuracy. Minimum realization 
theory was originally developed by Ho and 


Kc I man [137], using Markov parameters (pulse 
rf sponse functions). Questions regarding 
minimum realization from various types of 
Ir put-output data and generation of a minimum 
partial real izatlon were studied by Tether 
[144], Silverman [145], and Rossen and Lapldus 
[147], Rossen and Lapldus [148] successfully 
appl led Ho-Kal man and Tether methods to chemi- 
cal engineering systems. 

A common weakness of the preceding schemes Is 
that the effects of noise were not evaluated. 
Among follow-up developments along similar 
I nes, Kung [156] presented another algorithm 
In conjunction with the s I ngu I ar val ue decompcr- 
s tlon technique to treat the presence of the 
noise. Under the Interaction of structure and 
control disciplines, the Elgensystem Real Iza- 
t on Algorithm (ERA) [61] was developed by 
J iang and Pappa for modal parameter Identlflca- 
t on and model reduction for dynamic systems 
from test data. Based on a similar approach, a 
f eq ue ncy-doma I n ERA and a recursive ERA also 
wore developed [168,170]. Thorough treatment 
O’ the effects of noise on the ERA- Identl f led 
mod a I parameters was presented by Juang and 
Poppa [167]. Correlation of several modal 
testing methods was derived by Juang [169] via 
system real Izatlon theory. 

Ladder or lattice filtering first appeared In 
the reviewed literature by Morf, et al. [155] 
as a recursive method for solving the linear 
I oast squares problem. The term " I adder" or 
"lattice" originated from the shape of the data 
flow diagram of the technique. Interpreted In 
tsrms of the correspond I ng state space real Iza- 
tlon, the ladder Implementation uses a state 
vector that has a diagonal covariance matrix. 
Recursive Identification using the ladder 
representation has been extensively studied by 
Lee and Morf [159] and Frfedlander [88]. 
Recentl y , the algorithms have been appl led to 
Identification of flexible structures by Sund- 
erarajan and Montgomery [162,164] and Wlberg 
[165]. 

CONCLUDING REMARKS 

The field of modal testing has expanded con- 
tinuously over Its 40-year history. This 
growth Is largely associated with corresponding 
Improvements In computer capabilities. These 
Increases In computer capabll Ity have permitted 
more accurate and complete testing and data 
analysis to occur. Algorithms and approaches 
-thought too extensive In the past are new 
feasible. Several new uses for modal test data 
have also evolved. In particular, experimental 
nodal data are now being used frequently to 
t 1 recti y solve the problem at hand, rather than 
cnly for refining a finite-element model, the 
-traditional use In the aerospace community, 
'hese direct uses of test data include predlct- 
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fng the effects of physical changes to the 
structure using an experimental modal model and 
developing hybrid anal yt Ical/exper Imental 
models of overall system dynamics. 

For active control of space structures, the 
experimental modal data can also be used di- 
rectly for control law design, once the final 
configuration of the system Is tested using the 
control sensors and actuators. This task will 
require complete and accurate I dentl f Icatlon of 
the system while In orbit and Is motivating 
further technology Improvements to ensure 
success. Complex, built-up structures. In 
particular, still pose a sign! f leant chal lenge 
to the best ground-based methodology now avail- 
able. Success with large space structures will 
demand the combined efforts of the control and 
structural dynamics disciplines. The solid 
theoretical and methodological foundations from 
the control field should be combined with the 
extensive experimental knowledge from the modal 
testing field. Additional work Is needed to 
better understand and correlate current tech- 
niques from both fields. A principal goal Is 
to find a common basis to explain and to select 
from the myriad of possible techniques. 
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